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Abstract. We show how the measured particle ratios at RHIC can be used to provide 
non-trivial information about the critical temperature of the QCD phase transition. 
This is obtained by including the effects of highly massive Hagedorn resonances on 
statistical models, which are used to describe hadronic yields. Hagedorn states are 
relevant close to T c and have been shown to decrease rj/s to the KSS limit and allow for 
quick chemical equilibrium times in dynamical calculations of hadrons. The inclusion of 
Hagedorn states creates a dependence of the thermal fits on the Hagedorn temperature, 
Tji, which is assumed to be equal to T Cl and leads to an overall improvement of thermal 
fits. We find that for Au+Au collisions at RHIC at ^/snn = 200 GeV the best square fit 
measure, % 2 , occurs at T c ~ 176 MeV and produces a chemical freeze-out temperature 
of 170.4 MeV and a baryon chemical potential of 27.8 MeV. 



1. Introduction 

Hagedorn states, heavy resonances that follow an exponential mass spectrum, have 
been recently shown to play an important role close to the QCD critical temperature, 
T c . They have been used to find a low rj/s in the hadron gas phase pQ, which nears 
the string theory bound rj/s = l/(4n) [2]. Calculations of the trace anomaly including 
Hagedorn states also fit recent lattice results well and correctly describe the minimum 
c 2 . near the phase transition found on the lattice pQ. Furthermore, estimates for the bulk 
viscosity including Hagedorn states in the hadron gas phase indicate that (/s increases 
near T c pQ. 

Recent calculations using Hagedorn states to drive hadrons into chemical 
equilibrium indicate that hadrons do not need to be born in equilibrium at RHIC 
[31 HI El E] . This comes about because Hagedorn states allow for very short chemical 
equilibration times, which were previously not seen using only binary collisions. 
Previously, it was assumed that anti-baryons and anti-hyperons were already in chemical 
equilibrium following the phase transition from QGP to the hadron gas phase because 
gluon fusion more efficiently produces strange quarks. However, reactions including 
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Hagedorn states may also reproduce anti-baryons and anti-hyperons within the lifetime 
of a fireball. 

Because of the importance of Hagedorn states close to T c it is possible that they 
would have an effect on thermal fits, which are computed within statistical models to 
reproduce hadron yield ratios in heavy ion collisions [7J [El El QUI EH E2J- For Au+Au 
collisions at RHIC at ^snn = 200 GeV, specifically, estimates for the chemical freeze- 
out temperature and baryon chemical potential range from T c h = 155 — 169 MeV and 
fib = 20 — 30 MeV [101 IH E2] . Moreover, because Hagedorn states are dependent on the 
limiting Hagedorn temperature Tjj = T c , a relationship between the chemical freeze-out 
temperature and the critical temperature can be found by including Hagedorn states in 
thermal fits [13]. This uniquely gives us the ability to distinguish between different 
critical temperature regions depending on the quality of the fit obtained using the 
statistical model. Thus, in this paper we explore the possibility of including Hagedorn 
states into thermal fits. 



In this paper we use a grand-canonical model to describe the particle densities from 
which we can calculate the corresponding ratios as described in detail in [8]. Hagedorn 
states are included in our hadron resonance gas model via the exponentially increasing 
density of states [H] 



where M = 2 GeV and mf. = 0.5 GeV. We consider two different lattice results 
for T c : T c = 196 MeV [T^l US] (the corresponding fit to the trace anomaly is then 
A = 0.5GeV 3/2 , M = 12 GeV, and B = (340MeK) 4 ), which uses an almost physical 
pion mass, and T c = 176 MeV [17] (the corresponding fit to the energy density leads to 
A = 0.1GeV 3/2 , M = 12 GeV, and B = (SOOMe^) 4 ). Both are shown and discussed in 
[6]. Furthermore, we take into account repulsive interactions using volume corrections 
[DEI HE] that are thermodynamically consistent. Note that B is a free parameter based 
upon the idea of the MIT bag constant. Also, we include all the known hadrons below 
with mass below 2 GeV in our analysis. Note that B is a free parameter based upon 
the idea of the MIT bag constant. Note that due to the volume corrections Tjj > T c for 
a discussion on this see [HI [19]. 

In our model we do not just consider the direct number of hadrons but also the 
indirect number that comes from other resonances. For example, for pions we consider 
also the contribution from resonances such as p's, a/s, and etc. The number of indirect 
hadrons can be calculated from the branching ratios for each individual species in the 
particle data book [20J . Moreover, there is also a contribution from the Hagedorn states 
to the total number of pions, kaons, and so on as described in [H [6]. Thus the total 
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number of "effective" pions can be described by 

N w = N n + J2 N i(ni) (2) 

i 

whereas the total number of "effective" protons, kaons, or lambdas (generalized as X) 
can be described by 

N x = N x + Y,N t .(X t ) (3) 

i 

where (X) is the average number of X = p's, K's, or A's. Here N is the total number 
of each species and (n$) is the average number of pions that each Hagedorn state 
decays into. To determine (X) we use the multiplicities in Fig. 2 of Ref. [3] from 
the microcanonical model in [21] such that 



p = 0.058 rrii - 0.10 
K + = 0.075 nii + 0.047 

A = 0.04 rrii - 0.07. (4) 

Clearly, they are all dependent on the mass of the i th Hagedorn state. 

In order to get an idea of the quality of the thermal fits, we define x 2 as 

( TD ex P jDtherm\ 2 

i * 

where Rj herm is our ratio of hadron yields calculated within our thermal model whereas 
R^ xp is the experimentally measured value of the hadron yield with its corresponding 
error o\. Then, and T c h are varried until the minimum x 2 is found. We use the 
experimental data from both STAR [22] and PHENIX [23] at mid-rapidity for Au+Au 
collisions at RHIC at v /s/vjv=200 GeV. Specifically, we observe the ratios: n~ p/p, 
K~/K + , K+/7T+, p/n + , and (A + A)/tt+. All of which are calculated by STAR [22]. 
However, only tt~/tt + , p/p, K~/K + , K + /tt + , p/ir + are given by PHENIX. Because 
there is such a significant difference between p/tt + from PHENIX and STAR we choose 
only the value from STAR so that we can compare are results to [UJ where they also 
exclude p/ir + from PHENIX. It should be noted that Ref. pj] includes more ratios than 
we do such as multi-strange particles and resonances, which are not included in this 
paper. This is because the purpose of this paper is not to confirm their results, which 
have already been confirmed in [12], but rather to compare thermal fits that include the 
contribution of Hagedorn states and those that exclude them. 



3. Results 

We show the thermal fits for a hadron gas excluding Hagedorn states in Fig. [TJ where 
T ch = 160.4 MeV, /x 6 = 22.9 MeV, and x 2 = 21-2. This is almost identical to [TT] where 
T ch = 160.5 and p b = 20 MeV. 

The inclusion of Hagedorn states is our primary interest. Starting with the fit for 
the RBC-Bielefeld collaboration, we obtain T c h = 165.9 MeV, /i& = 25.3 MeV, and 
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Figure 1. Thermal fits for Au+Au collisions at RHIC at ^snn = 200 GeV without 
Hagedorn states with mass above 2 GeV. 

X 2 = 20.9, which is shown in Fig. [2} The x 2 is actually slightly smaller than in Fig. 
[TJ When we consider the lattice results from BMW, which are at the lower end of the 
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Figure 2. Thermal fits including Hagedorn states with mass above 2 GeV for Au+Au 
collisions at RHIC at ^snn = 200 GeV. 

critical temperature spectrum where T c = 176 MeV, we find T ch = 172.6 MeV, ^ = 39.7 
MeV, and x 2 — 17.8. The lower critical temperature seems to have a significant impact 
on the thermal fit. The lower x 2 lS due to the larger contribution of Hagedorn states at 
at T c h = 172.6 MeV, which is much closer to T c . 

The difference in the x 2 ' s f° r BMW and RBC-Bielefeld collaboration is directly 
related to the contribution of Hagedorn states in the model. Because the RBC-Bielefeld 
critical temperature region is significantly higher than its corresponding chemical freeze- 
out temperature the contribution of the Hagedorn states is minimal ~ 4-11% whereas 
the contribution from the Hagedorn states is 30 — 40% for T# = 176 MeV (see [32]). We 
find that the inclusion of Hagedorn states should not only provide a better fit but they 
also affect the chemical freeze-out temperature and the baryonic chemical potential. 
The more mesonic Hagedorn states are present the larger ^ becomes. Furthermore, our 
fits also have higher T^'s than seen in the fit without the effects of Hagedorn states. 
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4. Conclusion 

Our results for thermal fits without Hagedorn states concur well with other thermal fit 
models [IT] where the chemical freeze-out temperature (T C h = 160.4 MeV) is almost 
identical and the baryonic chemical potential = 22.9 MeV) is only slightly larger. 
The thermal fit with the known particles in the particle data group provides a decent fit 
with x 2 = 21.2. However, the inclusion of Hagedorn states provides an even better 
fit to the experimental data. In fact, we find \ 2 = 17.8, T c h = 172.6 MeV, and 
fib = 39.7 MeV for the BMW collaboration while for the RBC-Bielefeld collaboration 
we obtained \ 2 — 20.9, T ch = 165.9 MeV, and fib = 20.9 MeV. This provides further 
evidence [JJ HJ |5j E] that Hagedorn states should be included in a description of hadronic 
matter near T c . Since the chemical freeze-out temperature was found to increase from 
160 MeV to roughly 165 MeV (RBC-Bielefeld) or 172 MeV (BMW) when including 
Hagedorn states, this exemplifies the degree of uncertainty in extracting chemical freeze- 
out thermodynamical parameters by means of such thermal analyzes. 
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